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In ref. [Physica A 296, 405 (2001)], starting from the one parameter deformation of the expo- 
nential function exp^ ^{x) — + k?x'^ + kx^ ^ a statistical mechanics has been constructed 
which reduces to the ordinary Boltzmann-Gibbs statistical mechanics as the deformation parameter 
K approaches to zero. The distribution / — exp^^^ [— (iE + (3^ obtained within this statisti- 
cal mechanics shows a power law tail and depends on the non specified parameter /3, containing 
all the information about the temperature of the system. On the other hand, the entropic form 
Sk = J d^p{cK.f^^"' + C-k/^""*), which after maximization produces the distribution / and re- 
duces to the standard Boltzmann-Shannon entropy So as k — » 0, contains the coefficient Ck whose 
expression involves, beside the Boltzmann constant, another non specified parameter a. 

In the present effort we show that Sk is the unique existing entropy obtained by a continuous 
deformation of So and preserving unaltered its fundamental properties of concavity, additivity and 
extensivity. These properties of Sk permit to determine unequivocally the values of the above men- 
tioned parameters /3 and a. Subsequently, we explain the origin of the deformation mechanism 
introduced by ac and show that this deformation emerges naturally within the Einstein special rela- 
tivity. Furthermore, we extend the theory in order to treat statistical systems in a time dependent 
and relativistic context. Then, we show that it is possible to determine in a self consistent scheme 
within the special relativity the values of the free parameter k which results to depend on the light 
speed c and reduces to zero as c — > oo recovering in this way the ordinary statistical mechanics and 
thermodynamics. The novel statistical mechanics here presented preserves unaltered the mathemat- 
ical and epistemological structure of the ordinary statistical mechanics and is suitable to describe 
a very large class of experimentally observed phenomena in low and high energy physics and in 
natural, economic and social sciences. Finally, in order to test the correctness and predictability 
of the theory, as working example we consider the cosmic rays spectrum, which spans 13 decades 
in energy and 33 decades in flux, finding a high quality agreement between our predictions and 
observed data. 
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I. INTRODUCTION 

The following one-parameter deformations of the ex- 
ponential and logarithm functions 



^^{X) - [^/l + K^ +KxY'^ , (1.1) 



exp 



2k 



(1.2) 



which reduce to the standard exponential and logarithm 
respectively as the real deformation parameter n ap- 
proaches to zero, has been introduced recently in Ref. [Q . 
The above functions have many very interesting proper- 
ties 0, |[ ^, ^] (some being identical to the ones of the un- 
deformed functions) that permit to construct a statistical 
mechanics (and thermodynamics) which generalizes the 
standard Boltzmann-Gibbs one. This K-deformed statis- 
tical mechanics preserves unaltered the structure of the 
ordinary one and can be used to explain a very large class 
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of experimentally observed phenomena described by dis- 
tribution functions exhibiting power law tails. The ar- 
eas where this formalism can be applied include among 
others, low and high energy physics, astrophysics, econo- 
physics, geology, biology, mathematics, information the- 
ory, linguistics, etc. §|0||. 

In ref. [1| it has been shown that the statistical distri- 
bution 



/ 



expj^j(-/3[ii;-/i] 



(1.3) 



which generalizes the Boltzmann-Gibbs distribution, can 
be obtained also by maximizing, after properly con- 
strained, the entropy 



5«= / d"z;(c,/i+'' + c_,/i-«) 



(1.4) 



which reduces to the standard entropy S'o as the defor- 
mation parameter approaches to zero. The coefficient c^, 
which also absorbs the Boltzmann constant , depends 
on a free parameter a (see Eq. (65) of ref. |jl|) which 
remains to be determined together with the parameter (3 
which contains the information about the temperature T 
of the system. 
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A first question which arises naturally is if it is possible 
and how to find any criterion which allows us to fix the 
parameters /3 and a or at least express these in terms of 
the deformation parameter n, in order to reduce the free 
parameters of the theory. 

A second question regards the properties of the en- 
tropy 5'k. It is well known that the Boltzmann-Shannon 
entropy 5*0 is concave, additive and extensive. We know 
that Sk, is concave with respect to the variable /, but 
what happens about its additivity and extensivity? More 
in general, beside the Boltzmann-Shannon entropy other 
concave, additive and extensive entropies do exist? 

A third question is related to the physical mechanism 
which originates the deformation introduced by the pa- 
rameter K. In other words, does a more fundamental the- 
ory exist where this deformation emerges, or is it simply 
a purely mathematical tool? 

A fourth question is if it is possible to extend the the- 
ory originally developed in the framework of a classical 
kinetics to treat statistical systems in the context of a 
relativistic kinetics. 

A fifth question regards the deformation parameter k. 
This parameter will continue to remain free or is it pos- 
sible to determine its value self consistently within the 
theory? 

The present paper deals with the above questions and 
its purpose is double. Firstly, we will show that S^, is 
the unique existing, concave, additive and extensive en- 
tropy, beside the Boltzmann-Shannon entropy. As we will 
see these properties of are sufficient to determine un- 
equivocally the values of parameters /3 and a. Secondly, 
we will show that the deformation introduced by k is 
a purely relativistic effect and then we will explain the 
deformation mechanism within the Einstein special rel- 
ativity. Then, we will formulate a relativistic K-kinetics 
and we will calculate the value of k. 

Finally, in order to test the predictability and correct- 
ness of theory here proposed we will consider two sets of 
experimental data. In particular we will analyze the cos- 
mic rays spectrum which spans 13 decades in energy and 
33 decades in particle flux that, as it is widely known, 
violates the Boltzmann-Gibbs statistics. As we will see, 
we have a high quality agreement between the theory and 
the observed data. 

The paper is organized as it follows. In Sect. II, gener- 
alizing the approach proposed in ref.s |jl|, ^, we introduce 
a class of one parameter deformed structures and study 
their mathematical properties. 

In Sect. Ill, within this context, starting from the 
Jaynes maximum entropy principle we consider the most 
general class of deformed statistical mechanics preserving 
the main features of the standard Boltzmann-Gibbs one. 

In Sect. IV, we show that the entropy introduced 
in ref. |0 is the only one existing beside the Boltzmann- 
Shannon entropy Sq which is simultaneously concave, ad- 
ditive and extensive. Then the statistical mechanics and 
thermodynamics based on can be viewed as a natural 
extension of the Boltzmann-Gibbs one, recovering this 



last as the deformation parameter k approaches to zero. 

In Sect. V, we consider the mean properties of k- 
exponential and K-logarithm which have a fundamental 
role in the formulation of the new statistical mechanics. 

In Sect. VI, we extend the formalism to a time de- 
pendent and relativistic context. In particular after in- 
troducing the relativistic K-kinetic evolution equation we 
study its stationary state and prove the H-theorem. 

In Sect. VII, we explain the origin of the K-deformation 
and show that it emerges naturally within the Einstein 
special relativity. 

In Sect. VIII, we propose an approach which permits 
to determine within the special relativity the value of the 
parameter k. 

In Sect. IX, we compare two sets of experimental data 
with the predictions of the present theory. 

Finally in Sect. X, some concluding remarks are re- 
ported. 



II. DEFORMED MATHEMATICS 

A. Generator of the Deformation 

Let g{x) be an arbitrary real function of the real vari- 
able X, that we call generator of the deformation, having 
the following properties 

i) g{x) e C7-(R); 

ii) 9{~x) = -g{x); 

iii) dg{x)/dx > 0; 

iv) g{±(X)) — ±co; 

v) g{x) ~ X, for x --^ 0. 

Starting from the generator g{x), we construct the real 
function x^^y of the real variable x and depending on the 
real parameter follows 



a^j^j = — arcsinh(ji(Ata;) . (2.1) 

Its properties descent directly from the ones of the gen- 
erator g{x): 

i) eC-(R); 

ii) = -^w'^ 

iii) dxj^j /dx > 0; 

iv) (±oo)(^j = ±oo; 

v) x,. , « X, for a; — > and then 0,. , =0: 

vi) x,^, MX, for K — > and then x,„, = x: 

' x'^i {U> ' 

vii) X , , — x,^ , . 

' {-K.} {K} 

Together with the function x,^, one can introduce the 
inverse function x*"', through (a;*'^*)(^j = (a^j^j)*''^ = x, 
which assumes the form 

x^'^} ^^"^(sinh Kx) . (2.2) 
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B. Deformed algebra 

K 

Proposition 1: The composition law defined through 

{^®v)i^} = ^{«} ' (2-3) 

which reduces to the ordinary sum as k — > 0, namely 



X (B y = X + y, IS a. deformed sum and the algebraic 

K 

structure (R, ©) forms an abelian group. 

Proof: Indeed from the definition of x^^j the following 

K 

properties of follow 

1) associativity property: {x (B y) ® z = x ® {y (B z); 

K K 

2) neutral element: x(BO — 0(Bx — x] 

K K 

3) opposite element: x (—x) — {~x) a: = 0; 

4) commutativity property: x (B y = y (B x. 

K 

Of course the K-difference indicated with is defined 
as X Q y = X (B (—y)- 

Proposition 2: The composition law (8) defined through 



Proof: These prope rties follow directly from the defi- 
nitions given by Eq.s ( |2.3| ) and (^). 

Proposition 5: The function a;^^j and its inverse x^'^^ 
obey the following scaling laws: 



with 



X {ft'} — z X^^j 



K = k/z . 



(2.8) 
(2.9) 

(2.10) 
(2.11) 



Proof: These laws follow from the definitions of the 
functions x^^j and x^'^^ . 

Proposition 6: The pseudodistributive law 

z-{x(By) = {z-x) %\z-y) , (2.12) 

K 

holds and then the structure (R, 0,-) is a pseudofield. 

Proof: Indeed by using the propositions 1 and 5 one 
obtains 



(2.4) 



which reduces to the ordinary product as k — > 0, namely 



X IS) y = xy, is a, deformed product and the algebraic 

structure (R — {O},®) forms an abelian group. 

Proof: Indeed from the definition of x^^j we have the 
following properties 

1) associative law: (a; ® y) z = a; (g) (y ® z); 

2) neutral element: x®I — I®x — x\ 

K K 

3) inverse element: x®x — x®x = I; 

4) commutative law: x ® y — y ® x\ 

being the neutral element / = 1'^''^ while the inverse 
element of a; is S = (l/a:j^j )'f'*^. 

Of course the K-division is defined as x j/ = x y. 

Proposition 3: The deformed sum and product 
obey the distributive law 



z-{x®v)^ z-((x0y)(,, 



z ■ \ — X , H — y , 

' Z {K } Z {K } 



z-\\{x'By')^^,^ 
z • ( ( i (x' 2/') 



x' y' 

(z • x) ©"^ (z • y) 



C. Deformed derivative 



z (x y) = (z x) (z y) 



(2.5) 



and then the algebraic structure (R, ©,0) forms an 
abelian field. 



Proof: This proposition follows from the (2_^) and 

K K, 

We remark that the field (R, 0,0) is isomorphic with 
the field (R, +, •). Moreover z • (x y) 7^ (z • x) © (z • y) 

K 

and then the structure (R, ©, •) it is not a field. 

Proposition 4-' The function x*"^* has the two following 
properties 



x***' y 



{ft} 



(xy 



\{ft} 



(2.6) 
(2.7) 



Consider the two algebraic structures {X, ffi, • ) and 
(y, +, • ) with X = R and y = R. Let us introduce the 

set of the functions J" ^ {f : X ^ Y} with J" C C°°{X). 

(2.13) 



The K-differential d,^.,x is defined as 

ift> 



d,^, X — lim X z . 



and results to be 



(2.14) 



We define the K-derivative for the functions of the set 
T through 



(2.15) 
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with x^z € X and f{x),f{z) € Y. We observe that 
the K-derivative, which reduces to the usual one as the 
deformation parameter k — > 0, can be written in the form 



It resuhs that 



dfjx) ^ dfjx) 



1 



dfixl 

dx^i^Jdx dx 



(2.16) 



from which clearly it appears that the K-derivative is gov- 
erned by the same rules of the ordinary one. 



D. Deformed exponential 

The K-exponential exp^^^ {x) £ J- is defined as eigen- 
state of the K-derivative 



d exp^^j {x) 



dx 



exp,^j(a:) 



and is given by 
It results that 

^^P{0} = exp(a;) , 
exp^_^j {x) = exp^^j (x) 



(2.17) 

(2.18) 

(2.19) 
(2.20) 



The K-exponential, just as the ordinary exponential, 
has the properties 



exp^^j(x)eC°°(R); 


(2.21) 




(2.22) 


expj^j(-oo) = 0+; 


(2.23) 


exp^,^(0) = l; 


(2.24) 




(2.25) 


exp^^j(x)expj^j(-x) = 1 . 


(2.26) 



Furthermore, the K-exponential has the two properties 



= exp , irx) 



(2.27) 



exp^^j {x) expj^j {y) = exp^^^ {x ® y) , (2.28) 

with r G R, and can be expressed in terms of the gener- 
ator g{x) as 



'^^P{«:} (^) = +5('«a;)^ + 9{k.x) 



i/« 



(2.29) 



E. Deformed logarithm 



The K-logarithm In^^^ (x) is defined as the inverse func- 
tion of the of K-exponential, namely In^-^j (exp^^^ x) = 



ln{«,}(a;) = (Inx) 



(2.30) 



In 



{-K} 



(2.31) 
(2.32) 



The K-logarithm, just as the ordinary logarithm, has 
the properties 



V,(x)eC7-(R+); 


(2.33) 


^ V,(x)>0; 


(2.34) 


ln,,,(0+) = -(X3; 


(2.35) 


in,K,(i) = 0; 


(2.36) 




(2.37) 


- -ln(^j(a;) . 


(2.38) 



Furthermore, the K-logarithm has the two properties 

ln^,,(x'-)-rln^,,j(x) , (2.39) 

(x y) - In^,) (x) © In^^j {y) , (2.40) 

with r G R, and can be expressed in terms of the gener- 
ator g{x) as 



1 



(2.41) 



Eq. ( ^ ) defines a very large class of deformed log- 
arithms varying the arbitrary function g{x). These de- 
formed logarithms can depend on many other parameter 
(through the generator g(x)) besides the parameter k. 
We recall briefly that in literature one can find other one 
]lO|, 11 or two parameter deformations of the ex- 



ponential and logarithm functions. Anyway in following 
we w ill co nsider the deformed logarithms defined through 
Eq. (2.41) and depending only on the parameter k. 



F. Deformed trigonometry 

We define the K-hyperbolic sine and cosine 
sinh(^j(a;) = ^ [exp^^^ (x) - exp^^^ (-x) 
coshj^j {x) = i expj^j (x) + expj^j (-x) 
starting from the K-Euler formula 

^^'^{K} = cosh^^j (x) ± sinhj^j (x) 

It is straightforward to introduce the K-hyperbolic 
trigonometry, which reduces to the ordinary one as k — > 
0. For instance, the formulas 



(2.42) 
(2.43) 

(2.44) 



cosh^^j (x) - sinh^^j (x) = 1 
sinh^^j (x) 



tanh(^j(x) = 
cothj^j (x) = 



cosh(^j(x) ' 
cosh^^^ (x) 
sinli(«:}(2;) 



(2.45) 
(2.46) 

(2.47) 



5 



still hold true. All the formulas of the ordinary hy- 
perbolic trigonometry still hold true after properly de- 
formed. The deformation of a given formula can be ob- 
tained starting from the corresponding undcformed for- 
mula, and then by making in the argument of the hyper- 
bolic trigonometric functions the substitutions x + y — > 

X (B y, and obviously nx — > a; © x... © a; (n times). For 
instance it results 

sinti^^} {x © y)4«inh^^j {xQy) = 2 sinh^^j {x) cosh^^j (y) , 

(2.48) 



tanh(„j [x) + tanh(^, {y) 



sinh^^j {x © y) 



cosh(^j(a;)cosh^^j(y) 

(2.49) 

and so on. 

The K-De Moivre formula involving hyperbolic trigono- 
metric functions having arguments of the type rx with 
r g R, assumes the form 

[cosh(^j(x)±sinh(^j(a;)]' = cosh^^^^j(ra;)±sinh(^/^j(rx) . 

(2.50) 

Also the formulas involving the derivatives of the hy- 
perbolic trigonometric function still hold, after properly 
deformed. For instance we have 



d sinh^^^(a;) 
d tanh^^j {x) 



dx 



cosh(^j(a;) , 



[cosh^^j(x)]2 



(2.51) 
(2.52) 



and so on. 

The K-cyclic trigonometry can be constructed analo- 
gously. The K-sine and K-cosine is defined as 

sin^K} (x) = -« sinh^^j (ix) , (2.53) 
coS{^} (a;) = cosh^^j (ix) . (2.54) 

We remark that it results: sin,^,(a;) — sin(x . ) and 
cos.^i (x) = cos(a;, . , ). 

G. Deformed inverse functions 

The K-inverse hyperbolic or cyclic trigonometric func- 
tions can be introduced starting from the corresponding 
direct functions just as in the case of the undeformed 
mathematics. It is trivial to verify that K-inverse func- 
tions are related to the K-logarithm by the usual formulas 
of standard mathematics. For instance we have 



arcsiuj^j (a;) = — i In^^j {ix ^ v/ 1 — x^ , (2.55) 

1 1 + a; 

arctanh(^j(a;) = -In^^j -— , (2.56) 



H. Deformed product and sum of functions 

Let us consider the set of the non negative real func- 
tions V = {/, h, w, ...}. 

Proposition 7: The composition law (g) defined through 



(2.57) 



which reduces to the ordinary product as k — > 0, namely 
f ^ h = f ■ h, is a deformed product and the algebraic 

structure (T> — {0}, ®) forms an abelian group. 

K 

Proof: Indeed this product has the following properties 

1) associative law: (/ (8i h) w — f ^ {h ^ w); 

2) neutral element: / l = l{g)/ = /; 

K, K 

3) inverse element: / (g) (1//) = (1//) 0/^1; 

4) commutative law: f ® h = h® f . 

K Hi 

Of course the division can be defined through /0 h — 

K K 

f (E){l/h). The deformed K-power f®^ is defined through 



r^ = exp,^,(r V,/) , 



(2.58) 



and generalizes the ordinary power /''. In particular, 
when r is integer one has = / /... /, (r times). 

K K 

Proposition 8: The algebraic structure (I?,0) forms 
an abelian monoid. 

Proof: Indeed the element does not admit an inverse 
element. 

Furthermore, just as in the case of the ordinary prod- 
uct, it results / 0O = O0/ = O. 

Proposition 9: The composition law © defined through 

ffh = exp^^j |ln exp (in,^, /) + exp (in^^, h) } , 

(2.59) 

which reduces to the ordinary sum as the deformation 
parameter approaches to zero, namely / ©/i = / + /i, isa 



deformed sum and the algebraic structure {V, 0) forms 

K 

an abelian monoid. 

Proof: Indeed this sum has the following properties 

1) associative law: {f (B h) (B w — f (B {h® w); 

2) neutral element: / ©0 = 0©/ = /; 

K K 

3) commutative law: f (B h = h (B f ■ 

K K 

We remark that the product and sum © are dis- 
tributive operations w (/ © /i) = (w /) © (w /i). 

The product allows us to write the following prop- 
erty of the K-exponential 

Equivalently Eq. ( ^.60 ) can be written also in the form 
ln^.}(/® M-ln^«j(/) + ln(,,(M . (2.61) 



and so on. 



Eq. (2.61) gives a relevant property for the K-logarithm. 
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Finally, starting from the definition of the K-power 
Z®'", we obtain the following relation 



(2.62) 



The relations given by Eq.s ( ]2.6lD and ( |2.62D , which 
express two mathematical properties of the K-logarithm, 
will be very useful in the next section in defining a new 
additive and extensive entropy. 



Then the distribution n can be viewed as maximizing the 
information content L- 



Ik = J d p Ji^{n) ; Jf,{n) = ^ J l^f^, (n/a) dn , (3.7) 
under the constraints 

d^pE n^U , (3.8) 



III. THE JAYNES MAXIMUM ENTROPY 
PRINCIPLE 

Let us consider the following non normalized statistical 
distribution involving the k— exponential 



/ = exp 



{K} 



We write the real non specified parameter P as 



\k„T 



(3.1) 



(3.2) 



being A a new real parameter, the Boltzmann constant 
and T the temperature of the system. 

In the following it will be useful to introduce the dis- 
tribution 



a exp 



{K} 



E- 



(3.3) 



which is related with / through n = a f, being a another 
new real parameter which will be determined together 
with A in the following. 

We recall that in the ordinary statistical mechanics the 
mean value of a given physical quantity A{p,n) depend- 
ing on the variable p and the distribution n — n{p), is 
defined as 



< A{p,n)> = 



J (PpA{p, n) n{p) 
J (Pp n{p) 



(3.4) 



Analogously, in the case where A = A{pi,p2,ni,n2) de- 
pends on two independent variables pi , p2 and on the two 
independent distribution functions ni{pi) and n2{p2), we 
have that the mean value is given by 



<A> = 



J d^pid^p2 A ni{pi)n2{p2) 
J d^pi d^p2 ni(pi) 71-2(^2) 



It is easy to verify that the stationary distribution n 
can be obtained as solution of the following variational 
equation 



J d^p -fcs A J In^^j (n/a) dn-^En+j;: 







p n 



N 



(3.9) 



imposing the conservation of the mean energy and of the 
particle number respectively. Note that the chemical po- 
tential fi should be chosen in such a way to set the particle 
number equal to unity namely = 1 

We observe that when k = it results Jo (n) — nln n 
and the information content Iq is the mean value of the 
ordinary logarithm. In this case the above variational 
equation expresses the Jaynes maximum entropy prin- 
ciple which conducts to the Boltzman-Gibbs statistical 
mechanics. In the following, in analogy with the stan- 
dard statistical mechanics, we require that must be 
expressed as the mean value of In^^jTi. To do so we must 
consider the subclass of the deformed logarithms obeying 
the condition 



A 



ln{^)('^/a) dn = n In^^^i 



(3.10) 



The above condition, as we will see in the following, al- 
lows the simultaneous determination of both the form of 
K-logarithm and the values of the free parameters a and 
A. Firstly we observe that this class contains the stan- 
dard logarithm Inn, for which results k = 0, a = 1/e, 
A = 1. It will be the task of the next section to inves tigat e 
on the existence of new additional solutions of Eq. ( ^.10 ) , 
beside the standard logarithm. Taking into account this 
condition we can write the variational equation ( |3.6| ) in 
the form 



S 

Sn 



J d^p (^-k^ n In^^jn - 



n+^n 1=0 



We define the K-entropy through 



(Pp n In,. 



{K} 



(3.11) 



(3.12) 



(3.5) so that 5k can be viewed as proportional to the mean 



value of the In^^^n, namely 



Sk. — k^ 



< ln,^.,n> 



(3.13) 



In this definition of Sk, we have a perfect analogy with 
the Shannon entropy which is the proportional to the 
mean value of the Inn. It is remarkable that in both the 
definitions of 5*^ and So appears the standard mean value 



(3.6) given by Eq. (3.4) 
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Eq. (3.11) assumes the form 



^(-fcB< V,«>-^ <i?>+^ ) =0, (3.14) 



and then 







(3.15) 



The above variational equation can be viewed as defining 
a maximum entropy principle analogous of the Jaynes 
principle of the standard Boltzmann-Gibbs statistical 
mechanics We remark that this maximum entropy 
principle, in the form given by Eq. (3.11), holds only 
and exclusively for the subclass of k- logar ithms, which 
are solutions of the integral equation ( p.lO| ). 

We sh ow n ow that the families of entropies, defined 
through ( 3.12| ) and involving the K-logarithms which are 
solutions of Eq. ( |3.10D , have two important properties 
typical of the S hann on ent ropy. To do so we consider 
the properties ( 2.61 ) and ( 2.62 ) of K-logarithm, which 
rearrange as 



In 



r In^^jn 



(3.16) 
(3.17) 



with ni2 — ni (8)71,2 and n* — n^"^ . When the systems 

K 

1 and 2 described through ni and n2 respectively are 
statistically independent, and afte r ta ki ng i nto account 
the definiti ons o f the mean values (3.4), (3.£) and of the 
K-entropy ( 3.13 ), the two above properties of K-logarithm 
transform into the following properties for the K-entropy 



(3.18) 
(3.19) 



Eq.s ( |3l8| ) and ( p9| ) say that the entropies S^. defined 
starting from the K-logarithms which are solutions of 
Eq. ( |3.10 ) are additive and extensive just as the Shan- 
non entropy. The distribution ni2 describes the compos- 
ite system obtained starting from the systems 1 and 2 
while n* describes the scaled system related to the sys- 
tem described through n. Note that the state described 
through the distribution ni2 is different with respect to 
the state described through the distribution 7iirt2 result- 
ing S'«,[7ii2] < Si^[nin2]. 

Finally, from the concavity property of the deformed 
K-logarithm the concavity of S'^, follows 

5«[<ni + (l-i)n2] >t5«[ni] + (l-t)5«H ,(3.20) 

with < f < 1. 

As we have already noted, the ordinary logarithm is 
solution of Eq. ( 3.10 ) and then the Shannon entropy 



5*0 [n] = —kg f d^p n Inn 



(3.21) 



which is additive and extensive (ni2 = nin2 and n* 
n''), is admitted within the present formalism. 



In the next section we will show that Eq. ( [3.1C ) admits 
a new (only one) less evident solution. Then beside the 
Shannon entropy we have a new concave, additive and 
extensive entropy which is the K-entropy proposed in rcf. 

0- 



IV. THE NEW ADDITIVE AND EXTENSIVE 
ENTROPY 



We consider Eq. ( 3. 10] ) which, after performing a 
derivation with respect to n, assumes the form 



d 
In 



+ ln{^} n - A ln^^,(n/ a) = 



(4.1) 



In the following we will determine the explicit form of 
lUj^jTi by solving this differential- functional equation. 
We recall that In^^jn can be expressed in terms of the 
gene rator function according to Eq. (2.41) so that Eq. 
(4.1) becomes 



.9 

an 



2 / V 2 

(n/a)'" - (n/a)"" 



.(4.2) 



In the above equation the function to be determined is 
now the generator function g. To do so we make the 
following changes of variables 



i = K In n , 

z{t) — ^^"'^(sinht) , 

c — —Kin a , 



(4.3) 
(4.4) 
(4.5) 



so that Eq. ( [4.2| ) assumes the following simple form 

Kz'{t) + z{t)-Xz{t + c} = . (4.6) 

The property (v) of the generator g(x) imposes that 
z(t) obeys the two conditions z(0) — and z'{0) — 1. 
These conditions, if combined with Eq. (^^), can be 
equivalently written in the form 



Xz{c) = K 
Az'(c) = 1 



(4.7) 
(4.8) 



It is more convenient to take into account these condi- 
tions and write Eq. ( |4.6| ) under the form 



z{t + c) z{t) z'{c) + z'{t) z{c) 



(4.9) 



After recalling the prope rty (ii) of g{x), which imposes 
that z{—t) = —z{t), Eq. ( |4.9| ) can be written as 



2 z{c) z'{t) = z{c + t) + z{c - t) 



(4.10) 



Let us introduce the new function w(t) = z'{t). We can 
see that Eq. (4.10), after deriving with respect to c, 
transforms into the following functional equation 



2 w(c) w{t) = w(c + t)+ w{c - t) 



(4.11) 



8 



Finally the nonlinear transformation defined through 
w(t) = cosh^(t) permits us to write the last equation 
under the form 



cosh[c(c)+ew] + cosh [e(c) - ew] 

= cosh [^(c + t)] + cosh [^(c - t)] . 



(4.12) 



I t is t rivial to verify that the most general solution of Eq. 
(4.12) is given by 



m^rt , (4.13) 

with r an arbitrary real parameter. 

Shannon solution: We note that in the case r = we 
obtain g{x) = sinhx and then In^^j(n) = Inn. This is 
t he w ell known standard logarithm which, inserted in Eq. 
(3.12), produces the Shannon entropy. 

The new solution: We consider now the case r ^ 0. 
It is easy to realize that in this case the generator g{x) 
assumes the form 



g{x) — sinh 



— arcsinh (rx) 
r 



(4.14) 



For simplicity of the exposition we firstly discuss the case 
r — 1 for which g(x) = x. After some simple calculations 
we obtain that only when 



-1 < K < 1 



(4.15) 



it is possible to determine the real constants A and a 
obtaining 



1 - n 

1 + K 



1/2k 



(4.16) 
(4.17) 



In this case the generator g(x) — x imposes the following 
expressions for the deformed logarithm and exponential 



exp 



X — X 



2k 



(4.18) 

(x) = (Vl + ^2x2 + Kx) ^^'^ . (4.19) 



For the general case where r ^ 1 we obtain the same so- 
lution given by Eq.s ( 4.15| )-(4.1E) with the only difference 
that in place of k now the scaled parameter r/t appears. 
Then we can set r = I without loosing the generality 
of the theory. Moreover for k = we obtain the Shan- 
non solution as a particular and limiting case of the new 
solution. 

Finally, after insertin g th e expression of the k- 

logarithm given by Eq. ( 4. IS ) into Eq. ( 3.12 ), we can 



write the new additive and extensive entropy in the fol- 
lowing simple form 



2k 



(4.20) 



We can write the entropy 5'^ also in terms of the dis- 
tribution / obtaining 

SM] = ~K J d'p {cj'+^ + c_./i-«) (4.21) 

where the coefficient c„ — a^^'^ /2k depends exclusively 
on the deformation parameter k and is given by 



1 

2fi- 



1 



(4.22) 



The above entropy is contained, as a particular case, in 
the class of entropies introduced previously in ref. 0] (In 
Eq. (65) of this reference it appears the non specified 
parameter a, while the Boltzmmann constant is absent 
because setted to have kg[l + k) a = 1). 



We recall that the entropy given by Eq. (4.20) is dif- 
ferent from the nonextensive entropy introduced in ref.s 
Jl^ , ITsl. Of course the statistical distribution defined 



through Eq.s (3.1) and (4.19), introduced previously in 
ref. ||l|, is also different from the distribution of the 
nonextensive statistics [|l^ and of the plasmas physics 



V. THE K-EXPONENTIAL AND k-LOGARITHM 

Let us report here the main mathematical properties, 
some of these reported in ref.s jl], ||, |^, of the func tions 
exp [x) and In^^j (x) defined through Eq.s (4.1?) and 



(4.19), respectively. 



We start by observing that the generator of t he d cfor- 
mation is the function g{x) = x and then from (2.1) and 
(p]2h we obtain 



^{K} 



— arcsinh [kx) , 

K 

1 



— sinh {kx) 



We have also (0, P- 58) 
1 1 3 



X F \ — , —: —: —K^x^ 
2 2 2 



K^X^ < 1 



(5.1) 
(5.2) 

(5.3) 

(5.4) 



which has a central role in quantum group theory ^ , 
results to be proportional to x^^^, namely we have the 
relation 



The following function 



which can be written also in the form 

ln{«}(e) 



(5.5) 



(5.6) 
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with q = e"" . We note that the well known symmetry 
of quantum group theory q <-> q~^ is related with the 
symmetry k <-^—k oi t he p resent the ory. We also observe 
that, exploiting Eq.s (|2.6|) and (2/7), we can obtain the 
two following properties of the function [ x ] 



[x + y]^ [x] e [y] , 
[xy] = [x] h [y] , 



(5.7) 
(5.8) 



with k' = {q- q^^)/2. 

Th e defi niti ons of K-sum and K-product given through 
(^.3[) and (2.4) respectively transform as follows 



X ® y 



— sinh f 



arcsinh (kx) + arcsinh (k?/) ^ , (5.9) 



K 1 . / 1 . 

X ® y = — sinh — arcsinh (kx) arcsinh {ny) 

K \ K 

' (5.10) 

In particular the K-sum assumes a very simple form 



X (B y = x\J 1 + K^j/^ + y\/ 1 + K^x^ 



(5.11) 



Starting directly from the K-sum given by (5.11), one 
obtains the following expressions for the K-diffcrcntial 
and K-derivative 



dx 



dx 



{K} 



dfjx) 
dx 



dfix) 



(5.12) 
(5.13) 



We consider now the functions exp^^^ (x) and In^^j (x) 
which can be written also as 



exp^^j X — exp ^— arcsinh Kx 
In^^j (x) = — sinh (k In x) . 

We remark the following concavity properties 
d^ 



dx 
dx 



2 exp(^j(x) > ; X e 



R 



-j^ ln,,,(x) < ; X >0 



(5.14) 
(5.15) 

(5.16) 
(5.17) 



A very interesting property of these functions is their 
power law asymptotic behavior 



expj^j {x) ^ I 2kx 



±1/|k| 



a;— >±oo 



1 



x^a+ 2|k| 
1 

^ 771 — r ^ 

x^+oc 2 \k\ 



(5.18) 
(5.19) 
(5.20) 



The Taylor expansion of K-exponential is given by 

X 

a„(,K) 

ri=0 



exp^^j(x) = ^a„(K) 



n 



K^x^ < 1 



(5.21) 



([0, p. 26), with the coefficients a„ defined as 
ao(K) = 1 ; ai{K) = 1 , 



J=0 



(5.22) 



It results that a„(0) = 1 and a„(— k) = a„(K). We note 
that the first three terms in the above Taylor expansion 
are the same as the ordinary exponential, namely 



exP(^j(x) 



1 



(5.23) 



In fig.l, the function cxjp ^^^{—x) for a fixed value of k 
is plotted. We note that the bulk of this function is very 
close to the standard exponential. Indeed the Taylor ex- 
pansion of exp^^j {~x) is the same, up to second order, 
of the one of cxp(— x). The tail of exp^^^(— x) behaves 
as a power law. Between the bulk and the tail an inter- 
mediate region whose extension depends on the value of 
K exists. 



10" 



10"' 



^ "lO"' • exponential 



O. 
<D 



10"^! 




region 



K= 0.3 



— exp (-x) 

exp(^j(-x) 



power law - 
region 



10 



10" 



10 



10 



FIG. 1: Plot of the function exp (—x) versus x for k = 0.3. 

{t} ^ ' 

This function is compared with the ordinary exponential and 
with a pure power law. 

In fig. 2, the function expj^j(— x) for some different 
values of k is plotted. We note that when k ^ the 
K-exponential approaches the ordinary one. The Taylor 
expansion of lnj^j(l + x) converges if —1 < x < 1 and 
assumes the form 



V}(i+^) = E^»(^) (-1)""' 



X" 

n 



(5.24) 



10 



10' 



Iff 



5? 10 



I 10^ 



10 



10" 





K = 1 .00 




K = 0.75 




K = 0.50 




K = 0.25 




K = 0.00 






10"^ 


10-^ 




FIG. 2: Plot of the function exp (—a;) versus x for some 
different values of k. The case k — corresponds to the 
ordinary exponential. 



(0, p. 25), with = 1 while for n > 1, 6„(k) is 

given by 

It results 6„(0) = 1 and bn{—K) = 6„(k). The first terms 
of the expansion are 



2 13 



(5.26) 



Another expansion involving the K-exponential (x^ < 
1) is the following 



n=0 



exp^^j {x) = exp ^ dn 
(0, p. 58), being 
d„ = 



(-1)" (2n)! 



(2n+l)22"(n!)2 



(5.27) 



(5.28) 



Exploiting this expansion, we can write the K-exponential 
as an infinite product of standard exponentials 

oo 

exp^^j(x) = []exp(d„/t2n^2„+i^ _ (5 29) 

On the other hand the K-exponential can be viewed as a 
continuous linear combination of an infinity of standard 
exponentials. Namely for Res > it results (0, p. 
1108) 

f°° 1 /x\ 
exP{„j i-s) = — J,/^ [-) exp i-sx) dx . (5.30) 



The following two integrals can be useful 

x'^~^ exp^^j(— x) dx 



[1 + (r - 



2)|«|] |2^|-'^ r(pt^-§) 



[1 1^1]^ _ ^2 r(^ + 



r(r) (5.31) 



1 / ^y^l 
ln,K,-j dx 



l + (r- 



DN r(^ 



r-l 
k| ' 2 



r(r). (5.32) 



We conclude the present section focusing our atten- 
tion to another interesting property of K-logarithm. We 
consider the following eigenvector equation 



V,{x)L{x)^l,L{x) , 



(5.33) 



and examine the case of the eigenvector L{x) = n^. It 
is trivial to verify that when 'Dq{x) = d/dx we obtain 
that the eigenvalue of this operator is the standard 
logarithm Iq = Inn. We pose now the question if it is 
possible to determine the operator V^i^x) associated to 
the same eigenvector L(x) = and having as eigenvalue 
the K-logarithm, namely 



= In^^, n 



(5.34) 



We obtain that this operator is the finite difference op- 
eration 



V,{x)L{x) 



L{x + k) — L{x — k) 
2k 



(5.35) 



which reduces to the standard derivative as the increment 
2k of the independent variable approaches to zero. 

One can find many other elegant and useful mathemat- 
ical properties for the K-functions which obviously we can 
not report here. 



VI. RELATIVISTIC KINETICS 

In the present section we treat the statistical system, 
considered previously in stationary conditions, within a 
relativistic and kinetic framework. The new relativistic 
kinetics here presented, in the limit c ^ oo, reduces to 
the classical kinetics considered in [|l[ . 

By using the standard notation of the relativistic the- 
ory we denote with x = x'^ = {ct, x) the four-vector 
position and with p ~ p" = (p^,p) the four- vector mo- 
mentum, being p^ = ^/p"^ + rrfic^ and employ the metric 
g^^- = diag {I, m. 
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Let us consider the following relativistic kinetic equa- 
tion 

op J p' Pi Pi 

a(/'®/0-«(/®/i)l , (6.1) 

where the distribution / is a function of the four-vectors 
X and p, nam ely / = f{x,p). We note that the left hand 
side of Eq. ( |6.l[) is the same of the standard relativis- 
tic Boltzmann equation but the collision integral in the 
right hand side results to be more complicated, contain- 
ing the deformed product ® and the arbitrary function 

K 

a{f) which we suppose to be positive and increasing. The 
factor G is the transition rate which depends only on the 
nature of the two body particle interaction. 

The above equation in the case k = and a{f) = / re- 
duces to the already known relativistic Boltzmann equa- 
tion describing the standard relativistic kinetics [20|. 
Clearly in the case k 7^ and a{f) 7^ /, the above 
equation describes a new relativistic kinetics, radically 
different from the standard one. 

We anticipate that this new relativistic kinetics, which 
we will consider here, defines a statistics resulting to be 
independent on the particular form of the function a(/). 

Steady states: We consider now the steady states of 
Eq. (6.1) for which the collision integral becomes equal 
to zero. Then we have 



(6.2) 



and after taking into account the property (2.61) of the 
K-logarithm, we obtain 



ln,.,/ + ln,,,/i=ln,,,/' + ln,,^/( . 



(6.3) 



This last equation represents a conservation law and then 
we can conclude that ln,„, f is a summational invariant; 
in the most general case it is a linear combination of 
the microscopic relativistic invariants, namely a constant 
and the four-vector momentum. In ref. pO| it is shown 
that in presence of external electromagnetic fields the 
more general microscopic relativistic invariant has the 
form {p" + qA^/c) constant, being Ui, the hydrody- 
namic four- vector velocity with V^Ui, = (? . Then we can 
pose 



Km / 



(p" + gA7c) Uy - mc^ - /i 



(6.4) 



Consequently we obtain the following stationary distri- 
bution 



/ = exp 



{K} 



(6.5) 



In the case k = this distribution reduces to the already 
known relativistic distribution [EOl . 



The above equilibrium distribution, in the global rest 
frame where — (c, 0, 0, 0) and in absence of external 
forces {A^ = 0), simplifies as 



/ = exp 



{K} 



Xk^T 



(6.6) 



and assumes the same form of the distribution (3.1). 

We remark that for E — jj. ^ Xk^T this distribution 
presents a power law behavior, namely 



/ 



1/k 



(6.7) 



being E^ = k^T y/l - /2k. 

In order to introduce explicitly th e de pendence on the 
velocity variable in the distribution ( |6.6| ), we consider the 
expression of the relativistic kinetic energy E = E{v) 



E = Vm2c4 



with p = mj{v)\' and the Lorcntz factor given by 

1 



7(z;) 



v/l-i;Vc2 



(6.8) 



(6.9) 



After defining i] — 11 /mc we write Eq. (|6.6|) as follows 



Xk^T 



[7(w) - 1 - 7?] 



(6.10) 



Note that in the region v ^ c this distribution assumes 
the form 



1 2 



(6.11) 



which, after setting Xk^T — 1//3, coincides with the non 
relativistic statistical distribution, proposed in ref. 

H-theorem: In the standard relativistic kinetics it is 
well known from the H-theorem that the production of 
entropy is never negative and in equilibrium conditions 
there is no entropy production. In the following we will 
demonstrate the H-t heor em for the system governed by 
the kinetic equation (6.1). For simplicity of the notation, 
hereafter we omit the letter n in the symbol of K-entropy. 
We define the four- vector entropy S — = (5°, S), in 
terms of the distribution n — af, as follows 



5"^ 



-fc. 



d^P V 1 
-^p nln^^jn , 



(6.12) 



and note that S'"^ coincid es wi th the K-entropy defined 
previously through Eq. (4.20) while S is the entropy 
flow. If we take into account the relation ( [3.1C ), the 
above four-vector entropy can be written in terms of the 
distribution / as 



(6.13) 
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, d^p' (Ppi d?p\ <fp 

-K„AOL I -j: T\ -T\ 7- (_r 



It is trivial to verify that the entropy production d^S^ 
can be calculated starti ng f rom the definition of S'^ and 
the evolution equation (|6.l[), obtaining 



P'" Pi' P'l' 

-fc,Aam||j^(V^/)F^^. (6.14) 

Since the Lorcntz force F'^ has the properties p^Fi, = 
and dF'^ /dp^ = the last term in the above equation 
involving F'^ is equal to zero |M , namely 



d^p' d^pi d^p[ d^p 

7^ Pi° Pi" 7^ 



G 



a{r®f[)-a{f®h) ln,,j/. (6.15) 



Given the particular symmetry of the integral in Eq. 
(3.15) we can write the entropy production as follows 



- k„Xa 
4 



d^p' d^pi d'^p'i d^p 



G 



P'" Pi° P'i° P° 

' a{r®f[)-a{i®h) 

Finally, we set this equation in the form 



81,3" = - k„Xa 
4 



d^p' d^pi d^p\ d^p 



Pi" pf P" 



G 



K K 

ln,.j(/'®/0-ln,.>(/®/i 



(6.17) 



After imposing that a{h) is an increasing function, it re- 
sults [a[hi)-a{h2)] [ln^^j(/ii)-lnj^j(/i2)] > 0V/ii,/i2 and 
then we can conclude that 



d^S" > 



(6.18) 



This last relation is the local formulation of the rclativis- 
tic H-theorem which represents the second law of the 
thermodynamics for the system governed by the evolu- 
tion equation (6.1). 

Concerning the arbitrary positive and increasing func- 
tion a{f) appearing in the collision integral of the evolu- 
tion equation, we note that, if we suppose that obeys to 
the following condition 



a(/f/i)=a(/)a(/i) 
we recover the expression 

a{f) = exp(ln^^j /) 
proposed in ref. |^]. 



(6.19) 
(6.20) 



VII. PHYSICAL MEANING OF THE 
K-DEFORMATION 

In the present section we will show that the deforma- 
tion introduced by the parameter k emerges naturally 
within the Einstein's special relativity, so that one can 
see the K-deformation as a purely relativistic effect. 

Let us consider in the one-dimensional frame S two 
identical particles of rest mass m. We suppose that the 
first particle moves toward right with velocity vi while 
the second particle moves toward left with velocity 1)2 ■ 
The relativistic momenta of the particles are given by 
Pi —p{vi) and p2 =p{v2) respectively, being 



p{v) = 



v/1 - vyc^ 



(7.1) 



We consider now the same particles in a new frame S ' 
which moves at constant speed toward left with respect 
to the frame S. In this new frame the particles have 
velocities given by w 'j^ = v\ ®^W2 and w 2 = respectively, 
being 



1 + vw^jc^ 



(7.2) 



the well known relativistic additivity law for the veloci- 
ties. In the same frame S' the particle relativistic mo- 
menta are given by p 'j^ — p[v'-^ and P2 — 0, respectively. 
Up to now, we have simply recalled some well known 
concepts of the special relativity ||2l| . 

Let us pose the following question: if it is possible and 
how to obtain the value of the relativistic momentum p [ 
starting directly from the values of the momenta pi and 
P2 in the frame S. The answer to this apparently inno- 
cent question is affirmative. One, after straightforward 
calculations (see the theorem in this section), arrives to 
the following surprising result 



p[ ^ p{vi) <3) p{v2) ; K= — 



1 

mc 



(7.3) 



In words, the relativistic momentum p'l of the first par- 
ticle in the rest frame of the second particle is the n- 
deformed sum, with k = 1/mc, of the momenta pi and 
P2 of the particles in the frame S. 

Unexpectedly we discover that the K-sum is the addi- 
tivity law for the relativistic momenta. Eq.(7.3) which 
we write in the form 

K 1 

p{vi)®p{v2)^p{vi®''v2) ; K= — , (7.4) 

mc 

says that the K-deformed sum and the relativistic sum of 
the velocities are intimately related and reduce both, to 
the standard sum as the velocity c approaches to infinity. 
The deformations in both the cases are relativistic effects 
and are originated from the fact that c has a finite value. 
Eq.( [7.4[ ) follows as a particular case from the following 
theorem: 
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Theorem: Let 



Pi (Vt) 



(7.5) 



be the relativistic momenta of two particles (i = 1, 2) of 
rest mass mi and TO2 which move in the one dimensional 
frame S with speed vi and and V2 respectively. If we 
indicate with © the /c-sum defined through Eq. ( 5.11 ) and 
with (B'^ the velo city relativistic additivity law defined 
through Eq. (F^, it results that 



Pl (Ul) K P2 (W2) Pi{vi®''V2) 



1 



(7.6) 



mi 7712 fni c 

Proof: We start by using the definition of the K-sum, 
subsequently we use the explicit form of the relativistic 
momentum and finally we use the definition of the veloc- 
ity relativistic additivity law 



Pi (vi) If, P2 ("2) 




v/[i - ivi/cY] [1 
{vi (B" V2) 



- (^2/c)2] 
1 + VxV2l(? 



V/[1-(«1/C)2] [1 

Vi ®^ V2 




Pi {vi V2) 



(7.7) 



Trivially from Eq. (7^) one obtains Eq.([M) as partic- 
ular case, when mi = m2 = m. Note that the parameter 
K has different values in these two equations because the 
summed quantities in the two cases are different. 

We can easily explain the meaning of the deformed 
derivative. We indicate with p the relativistic momen- 
tum in the frame iS, and with dG/dp the derivative with 
respect to p of the Lorentz invariant scalar G. The same 
quantities in the frame S ' are indicated with p ' and 
dG/dp', respectively. It is trivial to verify that 

(7.8) 



dp 



{K} 



dp' 



and then we can conclude that the K-deformed derivative 
can be viewed as a standard derivative in an appropriate 
frame. 

In the next section we will consider, in the framework 
of the special relativity, the K-statistics of iV-identical 
particles, where k is a dimensionless parameter and we 
will determine its value. To do so, it is more convenient 
to write Eq. ([7.4|) in the form 



P («i) P (^2) _ P (vi V2) 



Kmc Kmc 
which holds for any k. 



(7.9) 



VIII. DETERMINATION OF THE 
PARAMETER k 

In the present section we calculate the value of the pa- 
rameter K which, due to the symmetry k — k of the 
theory, we consider positive, namely < k < 1. Clearly 
if we start from the K-statistical distributions /i and /2 
describing two independent statistical systems, we can 
construct the distribution /i /2 which describes a new 
composite system. This system in the case k = reduces 
to the one described through the distribution /i /2 . In 
the following we will assume that the distribution fi f2 
describes a state also in the case k ^ 0. Obviously this 
state is different from the one described by fi ® f2- As 
we will see, this simple but meaningful hypothesis is suffi- 
cient to determine the value of k in the case of relativistic 
systems. 

Taking into account the form of the distributions /i 



and /2 given by Eq. 
0), and the property 
write immediately 



(for simplicity we pose ^ — 
of the K-exponential one can 



El 



exp 



E^ 



Eo 



\k„T Xk„T 



(8.1) 



with Ei = Eivi). After some simple algebra we rewrite 
Eq. (B.l) as follows 



exp 



El 



exp, 
exp 



E. 



E^ 



(8.2) 



.3) 



We assume now that the K-exponential in the right hand 
side of the Eq. (8.2) has the same structure of the one 
given by Eq. (6.6). Clearly we must impose that Eq be 



with E3 = El + Ell El + E2 y/l + Ef/El and 

Xk^T 

En — . 
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exclusively expressed in terms of non statistical parame- 
ters. In the following w e w ill show that Eq — mc^. To 
do this we exploit Eq. ( [7.9[ ) obtained within the special 
relativity. Starting fr om t his equation and taking into 
account the property ( 2.28| ) of the K-exponential we ob- 
tain 



exp 



{K} 



p{vi) 



exp 



{K} 



P{V2) 



Kmc 



exp 



{K} 



p {vi V2) 



Kmc 



(8.4) 



Recall that we wish to calculate the parameter k which 
has a value that does not depend on the particle energy. 
Then, without loosing generality, we can consider Eq. 
(3.4) in the ultrarelativistic region (v c) where it re- 
sults p{v) ~ E{v)/c : 



exp 



E{vi) 



Kmc^ 



exp 



{K} 



exp 



E{vi < 



E{V2) 

Kmc^ 



Kmc'^ 



(8.5) 



After comparing Eq.s ( p.2[ ) and (3.5), one obtains the 
relation 



Km(? — \ k„T 



(8.6) 
only if 

Eq.(^.6D can be also written 



which is the same as the one given by Eq. 
we impose that Eq = m 
as 



•7) 



and results to be formally similar {k^T/c <-> p , k <-> 
v/c) to the relation defining the relativistic momentum 
given by Eq. (0). We can extract finally the value of k 
obtaining 



- = 1 



fc„T 



It is important to emphasize that this expression of the 
parameter k holds only under the above mentioned hy- 
pothesis and imposes that |k| < 1. This condition on 
t he v alue of k coincides with the one expressed by Eq. 
(4.15) and obtained in a completely different way. We 
have K — Q only if T = or if c = 00. The limiting case 
At = 1 is obtained if T = 00 or if mc^ = 0. 



At this point one can write the distribution (3.6) in 
the form 



/ = exp 



{K} 



1 E-^i 



K mc^ 



(8.9) 



Note that the statistical information of the system, 
namely the temperature is hidden exclusive ly i n the pa- 
rameter K. When E ^ 00 the distribution ( p.9| ) shows a 
power law asymptotic behavior 



^ ' '2E 



(8.10) 



The distribution ( |8.9| ) viewed as a function of the ve- 
locity becomes 



/ = exp 



7(1;) - 1 - ?/ 



Concerning its derivative one obtains 



V < c 



dv 

and then for 
dv 



r 



v/l + (7 - 1 - 



it results 



1 / V \ 



-1 + 1/2K 



(8.11) 



.12) 



.13) 



Then for k < 1/2 one has both f — and df /dv = Q in 

V — c. For K > 1/2 results / = and df /dv — ~oo in 

V = c. Finally for k ~ 1/2 results / = and df /dv = 
— l/2c in u = c. 

In the non relativistic region for which v c we have 



/ ~ exp^^j 

while in the limit c — * 
Maxwellian distribution 

Im = exp 



.14) 



Kmc^ 

00 one recovers the standard 



1 9 
■^mv 



k„T 



.15) 



The explicit form of the distribution (B.ll) when /i = 
—mc^ simplifies as follows 



/ = 



^l-i;2/c2 



i/« 



^l + y^2-vyc^ J 

and in the limit c 00 becomes 



V < c 



exp 



2fc„T 



(8.16) 



.17) 



In fig. 3, the distribution function given by Eq. (8.16) 
(after normalization) versus v/c for different values of k 
and then for different values of mc^/k^T according to 
Eq. 



, is plotted. 



IX. EXPERIMENTAL EVIDENCES 

Since long time it is known that the cosmic rays spec- 
trum, which extends over 13 decades in energy, from a 
few hundred of MeV (10* eV) to a few hundred of EeV 
(10^° eV) and spans 33 decades in particle flux, from 10"* 
to 10~29 [m"^ sr s GeV]~^ , is not exponential and then 
it violates the Boltzmann equilibrium statistical distri- 
bution cx exp(-£'/A;^T) ||, |||2|. Approximately this 
spectrum follows a power law E^"^ and the spectral index 
a is near 2.7 below 5 x 10^^ eV, near 3.1 above 5 x 10^^ 
eV and again near 2.7 above 3 x 10^* eV. On the other 
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FIG. 3: Plot o f the distribution function (after normalization) 
given by Eq. (S.16) versus vjc for some different values of k. 



hand it is known that the particles composing the cosmic 
rays are essentially the normal nuclei as in the standard 
cosmic abundances of matter. Then the cosmic rays can 
be viewed as an equivalent statistical system of identi- 
cal relativistic particles with masses near the mass of the 
proton (938 MeV). 

These above characteristics (relativistic particles with 
a very large extension both for their flux and energy) 
yield the cosmic rays spectrum an ideal physical system 
for a preliminary test of the correctness and predictability 
of the theory here proposed. 

We consider the statistical distribution /(i?) given by 
Eq. (|]|) or (U). The particle flux ^[E) cx v^f[E) 
can be calculated trivially if we take into account the 
relativistic expression linking E and p obtaining 



E 



1 



1 



(9.1) 

Note that this particle flux, in agreement with the obser- 
vational data, decays following the power law 



with 



(9.2) 



(9.3) 



Analogously the particle flux obtained starting from 
the Boltzmann-Gibbs statistical distribution is given by 



E 



exp 



E 



(9.4) 



We use these two theoretical distributions of particle 
flux to fit the cosmic rays data reported in ref. p3 . In fig. 



4, we show the observed data together with the theoreti- 
cal curves $(£■) (solid hue) and ^q{E) (dotted line). The 
curve $(£') corresponding io A — 10'"' [m? sr s GeV]~^ , 
mc^ = 938 MeV, = -375 MeV, and k = 0.2165 pro- 
vides a high quality agreement with the observed data. 
This agreement over so many decades is quite remark- 
able. From the value of k and mc^ and adopting Eq. 
~~7\} we obtain that k„T = 208 MeV. 
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FIG. 4: Plot of the cosmic rays flux versus energy. The solid 
line is the curv e ob tained within the present theory and is 
given by Eq. (|1|) with A = lO'^ [m^ sr s GeV]~^ , mc^ = 
938 MeV, fi = -375 MeV, k = 0.2165 and k^T = 208 MeF. 
The dotted line represents the theoretical curve obtained 
within the standard Boltzman-Gibbs statistics given by Eq. 
(0) with Ao = 1.3 • 10* [m^ srsGeV]-^, mc^ = 938 MeV 
and kgT — 208 MeV^. The observational data are collected 
by S.P. Swordy pi. 



In the same figure the curve ^o{E) (with Aq — 1.3 • 10* 



' sr s GeVy 



938 MeV and k^T ^ 208 MeV) 



which decays exponentially and can not fit the observed 
data violating the Boltzmann-Gibbs statistics, is re- 
ported . 

A short remark must be made at this point. Clearly, 
the power law asymptotic behavior of the spectrum ^(E) 
is imposed by the K-exponential whose origin is the k- 
sum. But the K-sum emerges naturally within the spe- 
cial relativity as the composition law of the relativistic 
momenta. Then we can conclude that the power law 
asymptotic behavior of the cosmic rays flux is simply the 
signature of the particle relativistic nature. 

It is widely known that the Boltzmann-Gibbs distri- 
bution exp{—E/kgT) originally proposed to describe a 
classical particle gas in thermodynamic equilibrium can 
be adopted to describe an enormous amount of phenom- 
ena in nature. On the other hand the power law tails 
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have been observed experimentally in several fields of sci- 
ence. Some times in particular fields, this power law has 
a name (e.g. Pareto law in econophysics, Gutenberg- 
Richter law in seismology etc). Furthermore the power 
law tail is preceded by an exponential region and be- 
tween the two regions exists a third intermediate region. 
It is worth remarking that the K-exponential defines a 
distribution which can describe simultaneously the three 
above regions (see fig.l) and then is particularly suitable 
to describe the above mentioned phenomena. 

As a working example we analyze the experimental 
data reported in ref. Q related to the rain events in 
meteorology. In fig. 5 is plotted the number density 
A'^ [events/ (year mm)] of rain events versus the event size 
M[mm] on a double logarithmic scale. We note that the 
data have a large extension (the abscissa spans 7 decades 
and the ordinate 5) and remark that its behavior is typi- 
cal of a class of experimental data which we find in several 
areas of science. In order to fit the experimental data we 
adopt the distribution 

7V = Aexpj^j(-/3M) , (9.5) 

and, as one can see in fig. 5, a remarkable agreement is 
obtained with A = 8 ■ 10* [events/ (year mm)], (3 — 75 
[mm^"'^] and k = 0.7. 



I 1 M M nij — r'-v-rTWVs| — ( rrTTWP«| — r~T~TTm^ — rnmrmt^ — t i a Ht t nj 




10' 10' 10' 10' 



Event size M [mm] 

FIG. 5: Plot of the number density [events/ (year mm)] of 
rain events versus the event size M[mm]. The solid line is the 
curve obtained within the present theory and is given by Eq. 
(U) with = 8 ■ 10* [events/(year mm)], /? = 75 [mm'^] 
and K = 0.7. The dotted line corresponds to the ordinary 
exponential function (k = 0). The experimental data are 
from ref. 

Clearly, one can hunt other mechanisms, different from 
the relativistic one, leading to K-statistics. Beside the im- 
portant problem of the agreement between the theoreti- 



cal curve and the observational data one does not neglect 
the epistemological problem concerning the structure of 
the theory which must be able both to explain the origin 
and to determine the value of any parameter appearing 
within the theory. 



X. CONCLUSIONS 

We summarize briefly the results obtained in the 
present effort. We have shown that beside the 
Boltzmann-Shannon entropy, the quantity 

i 

with In^^jX = {x'^ — x^'^)/2k and — 1 < k < 1, is the 
only existing entropy, simultaneously concave, additive 
and extensive. Starting from this entropy it is possible to 
construct a generalized statistical mechanics (and ther- 
modynamics) having the same mathematical and episte- 
mological structure of the Boltzmann-Gibbs one, which is 
recovered when the deformation parameter k approaches 
to zero. Within this generalized statistics the distribu- 
tion function assumes the form 

= a expj^j (-/?(£;,- ^) ) , 

with exp^^j(x) = (\/l + K^x^ + kx^ ^ while the con- 
stants a and /3 are given by a = [(1 — k)/(1 + k)]-'^/^", 
1/(3 = y/l — kj^T. The chemical potential fi can be 
fixed by the normalization condition. This distribution 
has a bulk very close to the exponential one while its tail 
decays following a power law rii cx -Ej ^^'^ . 

The origin of the deformation has its roots in the Ein- 
stein special relativity and the relativistic statistical sys- 
tems are governed by a kinetics obeying the H-theorem. 

We have shown that, within the special relativity, it 
is possible to determine the value of k, obtaining in this 
case 



so that the relativistic statistical meechanics does not 
contain free parameters. 

The theory can describe observational data in many 
fields. In particular we find a high quality agreement in 
analyzing the spectrum of the cosmic rays which violates 
manifestly the Boltzmann-Gibbs statistics. This is an 
important test for the theory because the cosmic rays 
are relativistic particles and their spectrum has a very 
large extension (13 decades in energy and 33 decades in 
flux). 
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